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Abstract 

We introduce three non-local observables for the two-dimensional Ising 
model. At criticality, conformal field theory may be used to obtain theoretical 
predictions for their behavior. These formulae are explicit enough to show that 
their asymptotics are described by highest weights hpg from the Kac table for 
c = ^ distinct from those of the three unitary representations (0, jq and ^). 

1 Introduction 

It is widely agreed that the Ising model in two dimensions is described at its critical 
temperature by a conformal field theory at c = |. The spectrum of the transfer 
matrix is that of the operator Lq © Lq and is described by a linear combination of 
the squared norm of characters of the three unitary representations at this value 
of the central charge, namely those with highest weights 0, ^ and |. The exact 
linear combination depends on the boundary conditions put on the geometry under 
consideration. 

These highest weights belong to Kac table hpg = (((m+1)]?— mg)^ — l)/4m(m+l) 
at c = 1 — 6/m(m + 1). For minimal models, the relevant (unitary) weights hpg are 
labeled by the integers p, q with l<p<m, l<g<m and 1 < p + q < 
m. As for the three unitary representations at c = |,m = 3 (corresponding to 
^00i^i.2 and /i2i), the Verma modules associated with the others hpg have singular 
vectors. The quotient by the subspace spanned by these vectors is irreducible but 
the inner product on the quotient space in not positive definite. Therefore these 
representations are not unitary and limits on field indices appearing in OPE are set 
to reject them. Prior to the work reported here, we did not know of any use for these 
non-unitary representations in the description of the Ising model, as for example in 
OPE's. 

We describe in each of the following sections an observable for the Ising model. 
They are somewhat unconventional as they are non-local objects. Using the tech- 
niques of conformal field theory, we are able to give predictions for their behavior 
and calculate their asymptotic behavior. The latter is described, in the three cases, 
by exponents hpg from Kac table that lead to non-unitary representations. 
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2 Crossing probability 
on Ising clusters 

The first observable is derived from one in percolation theory and we start by de- 
scribing it in this context. In percolation by sites on a square lattice, each site is 
declared open (closed) with probability p (resp. (1 — p)) independently of its neigh- 
bors. A configuration on a finite geometry is said to have a crossing between two 
disjoint intervals on the boundary if there is a cluster of open sites joining the two in- 
tervals. A common geometry is the rectangle, say of m x n sites, with the two disjoint 
intervals chosen to be the vertical sides. One then speaks of a horizontal crossing. 
A central quantity in percolation theory is the probability vr^'^'^™(r) of such crossings 
when the number of sites goes to infinity, the aspect ratio r = width/height being 
kept fixed. Cardy gave a prediction for this function vrJJ'^'^™(r) using conformal 
field theory and the agreement with numerical data is excellent 0. 

For the Ising model, we define 7r/i(r) = 7^]^^^^{r) as the probability of crossing on 
clusters of -|- spins, the limit on the number of sites being taken as above. Lapalme 
and one the authors have adapted Cardy's ideas to this case and obtained Monte- 
Carlo measurements to test their prediction Two steps in Cardy's reasoning 
cannot be extended straightforwardly to the Ising model and Lapalme and Saint- 
Aubin had to resort to one basic property of vr^'^'^'^" and vrjj'^™^, their scale invariance. 
If ttJj^™^ is described by a four-point correlation function, it must be that of a field 
of vanishing conformal weight. They therefore chose to use the second singular 
vector of height 6 in the Verma module with c = |, /i = 0. The ordinary differential 
equation obtained from this singular vector is of order 6 and the exponents are 0, | 
(twice degenerate), |, | and |. If one is willing to ignore the constraints on p and 
q, these exponents are precisely the conformal weights hpg of the fields appearing in 
the (naive) operator product expansion of 023 of conformal weight = 0. 

Lapalme and Saint- Aubin's argument is not as convincing as Cardy's which relies 
on better established ideas. The agreement of their prediction with numerical data 
is therefore welcome [Q. It is sufficient for the purpose of this letter to concentrate 
on the limiting behavior of 7ih{r) as r — >■ and r — »• oo. One expects 7r/i(r) — 1 
and — > in these two limits and the leading behavior prescribed by the ode is that 
of the exponent |. Because | is twice degenerate, a logarithmic behavior is allowed 
but appears to be ruled out by numerical data. The behaviors that are seen in the 
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Figure 1: Asymptotic behaviors of Hh as functions of r. The ordinates are log7r/i(r) 
for the top graph, log(l — 7r/i(r)) for the bottom one. 
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simulations are 

log7r/i(r) ^ ar + b, r — >• 

\og{l - TThil/r)) ^ - + d, r — > oo 
r 

with d = — 0.16647r and c = — 0.16657r. (The factor of vr stems from the change of 
variables between the aspect ratio r and the variable used in the ode. To obtain 
d and b, we used measurements of both horizontal and vertical crossings in ^ 
{TTv{r) = 1 — 7r/i(l/r)) and fitted the asymptotic behaviors above to the values of 
7ih for the 20 smallest (largest) aspect ratios r. See Figure ^) The measurement 
errors are of a few units on the fourth digit. There is little doubt that | = hs^s is 
the exponent describing these asymptotics. 



3 Contours intersecting the boundary of a cylin- 
der 

Consider now the description of the spin a at the boundary of a half-infinite cylin- 
der. Let 6i, 02, ■ ■ ■ , be the angles along the boundary where spinflips occur. Using 
Onsager's solution |^ or conformal field theory, it is possible to calculate the prob- 
ability density of configurations having precisely 2n flips. For example, this density 
s(^i, 62, O3, ^4) for the case n = 2 is proportional to 

(sin ^612 sin ^Osi)'^ - (sin ^613 sin ^024:)'^ 
+ (sin 16*14 sin^6'23)"^ 




Figure 2: Configurations contributing to different contour probabilities 
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if 61,62,9^, 64 appear in that order along the boundary and Oij is 6i — 6j. This den- 
sity does not reveal however which pairs {9i, 9j) are actually joined by the contours 
between same-spin clusters. For n = 1 there is only one possible pairing but, for 
n = 2, there are already 2. Using conformal invariance, we represent the half-infinite 
cyhnder by a disk (minus its center) and depict two possible distinct configurations 
(see Figure ^. Let l{9i, 62, O3, 64) (r(6'i, 62, 9^, 64)) be the density probability for the 
pairing of the left (resp. right) configuration. Two requirements allow for the deter- 
mination of / and r using conformal field theory. First their sum should reproduce the 
density when contours are ignored, namely l + r = s. Due to the singularity in s, it is 
natural to seek I within the solution space of the ordinary differential equation that 
describes the 4-point correlation function of the field 02,i of conformal weight |. Sec- 
ond, when 612 0, the probability Prob(^i^2|^3^4) = '/(^i, ^2, ^3, ^4)/s(^i, ^2, ^3, ^4) 
should go to 1 and Prob(6'26'3|6'46'i) = r/s to 0. 

These two requirements determine uniquely l,r and Prob(0i^2|^3^4) For the 
latter one gets 

Prob(0i02|^^3^4) = 

1 _ 9 ^(3) / .(5/3)/ N ^ f(5m 

2 2or(|)2/(°H-2) V 

where the anharmonic ratio is chosen to be 2; = (sin ^612 sin ^O^^) / (sin |6'i3 sin |6'24) 
and where f^\z) = l-z + z/{l - z) and f^^^\z) = zl 2i^i(-|, |, |, z)/{l - z). 
The behavior of this probability as 612 ^ (2; ^ 0) is 

ior(|)2 

9 r(i)^ +^[^ )■ 

The exponent | is the highest weight h^i of Kac table. 

The function Prob(6'i6'2|6'36'4) is plotted on Figure ^ together with Monte-Carlo 
measurements of this probability on a cylinder whose length is twice as long as its 
circumference. The four dots close to z = and the four close to z = 1 were measured 
on a cylinder with 32000 sites and samples were larger than 10^ configurations. 
Statistical errors on these 8 points are smaller than the size of the dots on the 
figure. (Details will be given in [Q].) 
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Figure 3: The function Prob(6'i6'2|6'36'4) as function of 2; G [0,1] with Monte-Carlo 
data. 

4 Homology class 

of Fortuin-Kasteleyn clusters 

The last observable describes the Ising model on a geometry without boundary. 
Denote by a (/?) a non-trivial cycle in the horizontal (vertical) direction on a torus 
of modulus T. Let a, 6 G Z be two integers with gcd(a, 6) = 1 and let 7r({a, 6}) be 
the probability that a spin configuration has a Fortuin-Kasteleyn cluster wrapping 
precisely a times around a and b around f3. (Configurations having simultaneously 
clusters of type {1,0} and {0, 1} are not included in the computation of n{{a,b}). 
These configurations are said to contain a cross.) In Figure ^the Fortuin-Kasteleyn 
clusters of three configurations on a square torus (r = i) are shown. The original 
signs of the Ising spins are depicted in black (— ) or in white (+). Only the leftmost 
configuration in this Figure contributes to 7r({l,0}). 

Using the Coulomb gas representation |^ it is possible to write an explicit 
formula for 7r({a, b}) (see [Q where this expression is given for various Potts models 
and where Monte- Carlo simulations checking it are reported). In particular for the 
Ising model 



where = Imr, q = e™'^ ^ the Qi are the elliptic theta functions and 77(g) is Dedekind 
function r]{(i) = Y['^=ii^ ~ 'f')- For Rer = and Xj — > 00, that is for a torus 



^r({(l,0}) 



1 g2(H-|g3(tr,)-g4(H| 



\rim I^2(r)| + |e3(r)| + |^4(r)| 
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Figure 4: Fortuin-Kasteleyn clusters for three configurations. Tliey contribute to 
7r{(l,0)}, 7r{(l, —1)} and 7r(cross) respectively. 



represented by a very narrow and tall rectangle, any configuration will almost surely 
contain a horizontal cluster and no vertical one. In that limit one expects that all 
7ir{{{a,b}) will vanish, including 7rT-(cross), except for 7rT-({l,0}). For that case 
q = e"'^'^' and one gets indeed as — > oo 

7r,({l, 0}) ^ 1 - (q')'sf,(q^) - {q^)\ f,{q^) 

where the /j are real analytic in their argument in a neighborhood of 0. The two 
leading terms are recognized to be twice the weights hi^2 and h^^^. (The doubling of 
highest weights in expansions in q^ is the natural thing to expect and accounts for 
the contributions of holomorphic and antiholomorphic sectors of the theory.) 

Of the three observables discussed in this note, this one (7rT-({l, 0})) is probably 
the less compelling. There is no differential equation here to dictate a finite set of 
exponents. The ". . . " in the above asymptotic expansion contains other exponents, 
namely other integral linear combinations of | and |. Not all these combinations 
however occur in Kac table. (We showed that the exponents arc restricted to the 
set {|, n = 0, 1, . . . , 7} U {| + |, n = 0, 1, . . . , 7} but we did not prove that all of 
these do occur.) 

Despite these comments, the fact that the highest weights hi2 and /133 describe 
the leading behavior seems remarkable. 



5 Conclusion 

What are the possible exponents of the Kac table at m = 3 if one identifies hpg and 
hp'q' whenever hpg — hp/g/ G Z? This amounts to asking the simple number theoretic 
question of which integers modulo 48 have a square root. The answer to the first 
question is the following list of /I's: the three unitary 0, and | and the non-unitary 



7 



|, |, ^, II and — ^. The question of whether other observables can be constructed 
that are ruled by ^, || and — ^ is more difficult. 

More to the point are the following questions. Which linear space of states 
must be considered so that non-local observables may be taken into account in the 
framework of conformal field theory? Why do statistical models seem to prefer 
unitary representations even when mutually non-local fields are considered (like the 
pair a-fi in the Ising model)? And why, when they do step out of the unitary 
representations, do they remain in Kac table? 



NOTE 

After the research reported here was completed, one of us (YSA) learned from 
Duplantier that various exponents for the Ising model have been introduced by him 
and his colleagues that do not belong to the small set {0, 3^^, |}. 

Here are two representative examples discussed in their work. For 0{n) models, 
consider the probability density of having L loops between two points x and y in 
the plane. Duplantier ^ shows that it decays as |x — ?/|~^^^ with xl = 2/1^/2,0- Not 
only do these exponents miss the unitary set {0, ^, |} but, for L odd, they are out 
of Kac table. This fact raises questions beyond the present letter. 

More recently Duplantier obtained the fractal dimensions of three proper- 
ties of Potts clusters, namely of the external perimeter [Dep), of the set of outer 
boundary sites (dimension Dh of the hull) and of the singly connecting sites that 
appear close under the scaling limit (Dsc)- For the Ising model he gets Dep = 
Dh = I and Dsc = -'^^ would be hard to believe that the | is a coincidence but 
the other two dimensions are out of Kac table. The role of these exponents for the 
argument presented here is unclear to us. 

The most intriguing connection with the present work has appeared recently. 
Read and Saleur [l^ consider nonlinear sigma models whose fields take values in 
supersymmetric coset spaces. They argue that, for the target space S"^"'^" (a su- 
persymmetric generalization of the sphere), the model shows a Ising-like transition. 
The spectrum of its conformal weights is described by their formula (3.7-8) with 
Co/ = I and gp = | and includes the set {0, ^, |} and the exponents observed 
here: | and |. The relationship between their models and the observables discussed 
here remains to be established. 
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